Simple solutions of relativistic hydrodynamics 
for cylindrically symmetric systems 



T. Csorgo/' 13 ' 1 F. Grassi, b ' 2 Y. Hama, b ' 3 T. Kodama c ' 

a MTA KFKI RMKI, H - 1525 Budapest 114, POBox 49, Hungary 
h IF-USP, CP. 66318, 05389-970 S&o Paulo, SP, Brazil 
C IF-UFRJ, CP. 68528, 21945-970 Rio de Janeiro, RJ, Brazil 



Abstract 

Simple, self-similar, analytic solutions of 1 + 3 dimensional relativistic hydrody- 
namics are presented for cylindrically symmetric fireballs corresponding to central 
collisions of heavy ions at relativistic bombarding energies. 
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1 Introduction 



The analytic resolution of 3 dimensional relativistic hydrodynamics is a dif- 
ficult task because the equations are non-linear and the transverse and the 
longitudinal motions of the fluid are coupled in a way that is rather compli- 
cated to handle. 



Historically, the first solutions of relativistic hydrodynamics were found for ide- 
alized, one dimensionally expanding sources: Landau's one- dimensional ana- 
lytical solution of relativistic hydrodynamics [1] and Bjorken's boost-invariant 
solution [2] , that correspond to ultrarelativistic bombarding energies. Both of 
these solutions are frequently utilized for various estimations of observables in 
ultra-relativistic nucleus-nucleus collisions. 
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The repertoir of relativistic hydrodynamics has been reviewed from the point 
of view of heavy ion physics in ref. [3]. Recently, Biro has found self-similar 
exact solutions of relativistic hydrodynamics for cylindrically expanding sys- 
tems [4,5]. Bjorken's one-dimensional solution can be naturally extended to 
a three-dimensional solution, assuming spherical symmetry instead of (1+1) 
dimensional longitudinal dynamics. Biro's solutions generalize the (1+3) di- 
mensional version of Bjorken's solution to the case of cylindrical symmetry, 
interpolating between the one-dimensional and the three-dimensional Bjorken 
solution. However, Biro's solutions are valid only when the pressure is inde- 
pendent of space and time, as happens at the softest point of the equation of 
state. 

The numerical solution of the 1+3 dimensional relativistic hydrodynamics is 
possible with presently available computers and methods, for example the 
variational method of smoothed particle hydrodynamics has been developed 
to study event-by-event fluctuations in relativistic heavy-ion collisions [6,7]. 

Here we describe a new family of exact analytic solutions of relativistic hy- 
drodynamics with a broad family of equations of state, assuming cylindrical 
symmetry of the particle emitting source. The physical motivation for this 
study is to consider the time evolution of central collisions in ultra-relativistic 
heavy-ion physics within the framework of an analytic approach. 

Our solution was motivated by the analytic solution of non-relativistic hy- 
drodynamics found by Zimanyi, Bondorf and Garpman (ZBG) in 1978 for 
low-energy heavy- ion collisions with spherical symmetry [8]. This solution has 
been generalized to the case of ellipsoidal symmetry by Zimanyi and collab- 
orators in ref. [9]. In ref. [10,11] a Gaussian parameterization has been in- 
troduced to describe the mass dependence of the effective temperature and 
the radius parameters of the two-particle Bose-Einstein correlation functions 
in high-energy heavy-ion collisions. Later it has been realized that this phe- 
nomeno logical parameterization of data corresponds to an exact, Gaussian 
solution of non-relativistic hydrodynamics with spherical symmetry [12]. The 
spherically symmetric self-similar solutions of non-relativistic hydrodynamics 
were generalized in ref. [13] in terms of an arbitrary scaling function for the 
temperature profile, and expressing the density distribution in terms of this 
scaling function. The spherically symmetric Gaussian solution has been gener- 
alized to ellipsoidally symmetric expansions [14], that provided simple analytic 
insight into the observables of non-central heavy-ion collisions [15]. The family 
of ellipsoidally symmetric solutions of non-relativistic hydrodynamics has been 
expressed in terms of a general scaling function for the temperature profile in 
ref. [16]. 

Our approach corresponds to a generalization of these recently obtained an- 
alytic solutions [12,13,15,16] of non-relativistic fireball hydrodynamics to the 
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case of relativistic longitudinal and transverse flows. In particular, an analytic 
approach, the Buda-Lund (BL) model has been developed to parameterize 
the single particle spectra and the two-particle Bose-Einstein correlations in 
high-energy heavy-ion physics in terms of hydrodynamically expanding, cylin- 
drically symmetric sources [17]. Here we attempt to find an exact family of 
solutions of relativistic hydrodynamics that may include the BL model as a 
particular limiting case. 

We attempt here to solve the full (1+3) dimensional relativistic hydrodynam- 
ical problem, in trying to overcome two shortcomings of Bjorken's well known 
solution. These two shortcomings of Bjorken's solutions are that i) it describes 
ultra-relativistic collisions, it is independent of scale parameter in the longitu- 
dinal direction and the rapidity distribution is flat; ii) it contains no transverse 
flow. 

We have described recently a new family of (1+1) dimensional solutions [18] 
that is able to describe an arbitrary rapidity distribution dn/dy by introducing 
scale-dependent quantities in the longitudinal direction, overcoming shortcom- 
ing i), but not addressing shortcoming ii). Here we consider the opposite case, 
overcoming shortcoming ii), the development of a scale dependent transverse 
flow in relativistic hydrodynamics, but not addressing shortcoming i). We ad- 
dress both shortcomings simultaneously in a subsequent manuscript [19]. 



2 The equations of relativistic hydrodynamics 



We solve the relativistic continuity and energy-momentum conservation equa- 
tion: 

d ll (nu»)=0, (1) 
d v T» v = . (2) 

Here the n = n(t, r) is the number density, the four- velocity is denoted by 
u n = U ^[t^ r ) — v ), normalized to u^u^ = 7 2 (1 — v 2 ) = 1, and the energy- 
momentum tensor is denoted by T^ v . We assume perfect fluid, 

r^f+pjuY-pf, (3) 

where e stands for the relativistic energy density and p denotes the pressure. 

We close this set of relativistic hydrodynamical equations with the equations 
of state. We assume a gas containing massive conserved quanta, 

e = mn + up, (4) 
p = nT. (5) 



3 



The equations of state have two free parameters, m and k. Non-relativistic 
hydrodynamics of ideal gases corresponds to the limiting case of m > T, 
v 2 < 1 and k — §. 

The energy-momentum conservation equations can be projected into a com- 
ponent parallel to w M and components orthogonal to u^, which are respectively 
the relativistic energy and Euler equations: 

u"V + (e + p)a/ = 0, (6) 
u v u^p + (e + p)u^u u -d u p = 0. (7) 



With the help of the equations of state and the continuity equation, the energy 
equation can be rewritten as an equation for the temperature, 

u^dfjT + -TdX = . (8) 

We adopt the following notational conventions: the coordinates are = 
(t, r) = (t,r x ,r y ,r z ), x^ = (t, —r x , —r y , —r z ) and the metric tensor is g^ u = 
= diag(l, -1,-1,-1). 

We solve 5 independent equations, the continuity, the (3 spatial components 
of) relativistic Euler, and the temperature equation, eqs. (1,7,8). The equa- 
tions of state, eq. (4,5) close this system of equations in terms of 5 variables, 
n, T and v = (v x ,v y ,v z ). 



3 Self-similar, cylindrically symmetric solutions 



As we are primarily interested in the effects of finite transverse size and the de- 
velopment of transverse flow, we assume that the longitudinal flow component 
is that of Bjorken's type, 

v z (t,r z ) = r -^. (9) 

A similar assumption has been made when a new family of (1+1) dimensional 
solutions of relativistic hydrodynamics was obtained in ref. [18]. However, in 
contrast to what has been done in [18], here we assume scale-invariance in the 
longitudinal direction. 

We search for self-similar solutions, that are scale dependent in the transverse 
directions, and depend only on the transverse radius variable r t = ^jr 2 x + r 2 

and the longitudinal proper time r z — Jt 2 — r 2 . Let us introduce the scaling 
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variable x as 



^2 I ^2 



* = ^, (io) 

and assume that, in the frame where v z = (longitudinal proper frame), 
the transverse motion corresponds to a Hubble type of self-similar transverse 
expansion, 

(12) 

where R = dR{r z )/dT z and hereafter we will designate by starred symbols 
the variables in the longitudinal proper frame. We assume that the scale R 
depends on time only trough the longitudinal proper time, t z . 

In a relativistic notation, the above form may be parametrized as 

■u M = (cosh ( cosh £, sinh£— , sinh£— , sinh(cosh£), (13) 

r t r t 

tanh£ = ||||r t =v; = -m, (14) 

cosh£ = - 1 . = 7t * , (15) 
'l - R 2 x 



coshC = — = 7; . (16) 
The space-time rapidity i] is defined as 

,.06Iog(|±i). (17) 
For a scaling longitudinal Bjorken flow we obtain 

C = V- (18) 

Using the above ansatz for the flow velocity distribution, we find that the 
continuity equation is solved by the form 



n(t, r x , r y , r z ) = n ( ^§ ) —^r-zG{x), (19) 



?2 x ^ 
\^ t z R 2 J cosh £ ' 



where G(x) is an arbitrary non- negative function of the scaling variable x 
and n , r z0 and R are normalization constants. We use the convention n = 
n(t , 0, 0, 0), r z0 = T z (t , r z0 ) and R = R(t z0 ), where r z0 is such that, together 
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with to, satisfies eq. (9). This implies that Q(x = 0) = 1. The temperature 
equation, eq. (8), is solved by 

T{t ,r M =T o {0^y K r(,). (20) 

The constants of normalization are T = T(t , 0,0,0) and JF(0) = 1. We find 
that the solution is independent of the form of the function J-{x), provided 
that F(x) > 0. 

Using the ansatz for the flow profile and the solution for the density and the 
temperature profile, the relativistic Euler equation reduces to a complicated 
non-linear equation that contains R, R and R and the variable x. Taking this 
equation at x = we express R as a function of R and R. Substituting this 
back to the Euler equation we obtain a transcendental equation for R 2 , and 
x. This equation has a particular solution if 

R = Rq. (21) 

In this case, the acceleration of the radius parameter vanishes, R — 0, and the 
solution is R = R + Rq(t z — t z0 ). Then the relativistic Euler equation reduces 
to 

l + £) (^ + 3R 2 ^J =2(l-xR 2 )[\ogg(x)F(x)}' , (22) 

where the lhs depends only on r z while the rhs is only a function of the variable 
x, hence both sides are constant. This implies that 

- = Ro, (23) 

thus Rq = RoT z o . Thus the origin of the time axis (fixed by the assumption 
of the scaling longitudinal Bjorken flow profile) coincides with the vanishing 
value of the transverse radius parameters. 

The solutions can be casted in a relatively simple form by introducing the 
proper time r, 



r = ^t 2 - r 2 = ^t 2 -r 2 x -r 2 y - r 2 z . (24) 

(25) 

Using this natural variable we find that 

v = (26) 

T fj. 

u» = —. (27) 
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Thus the velocity field of our solution corresponds to the flow field of the 
spherically symmetric scaling solution. However, in the scaling solution the 
temperature and the pressure distributions are dependent only on the proper 
time r, while in our case both the density and the temperature distributions 
are generally dependent on the scale variable x in the transverse direction. 

As the solution is relativistic, and it is defined in the positive light-cone, given 
by t > 0, we obtain a constraint for the transverse coordinate, r t < r z . This 
together with eq. (23), the solution for the scale R, implies that the scaling 
variable has to satisfy the constraint xR% < 1, which corresponds to the 
limitation that the velocity of the fluid can not exceed the speed of light. 

By replacing eqs. (21,23) into the Euler equation, eq.(22), one obtains 

log [(1 - xR 2 )^ 1+1 ^g(x)T(x)] = , (28) 

which gives, together with the condition ^(0)JF(0) = 1, 

Q(x)F(x) = (1 - i?^) -2(1+1/*). (29) 

In this family of solutions, the scaling functions for the temperature and the 
density distribution are thus not independent. However, a constraint is given 
for their product, hence one of them can be chosen as an arbitrary positive 
function. 

For clarity, let us introduce new forms of the scaling functions as 

T(x) =F(x) (1-R 2 x) 2 / K , (30) 
V{x) =g{x){l-R 2 x) 2 . (31) 

Then the constraint can be casted to the simple form of 

V(x)T(x) = 1. (32) 

This construction for the scaling functions of the transverse density and tem- 
perature profiles coincides with the method, that we developed for the solution 
of the relativisitic hydrodynamical equations in the (1+1) dimensional prob- 
lem, but here the transverse flow has a two-dimensional distribution, so the 
exponents and the scaling variables had to be re-defined accordingly. 

Let us summarize our new family of solutions of the 1+3 dimensional rel- 
ativistic hydrodynamics for cylindrically symmetric systems by substituting 
the results to the density, temperature and pressure profiles. 
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We obtain 



v = -, for Irl < t 
V 1 1 ~ 



(33) 




n T z 

n(t,r) = n (jyj V(x), 



Zl 



(34) 



(35) 



(36) 



(37) 



where po = UqT . Note that the scaling variable x is invariant for boosts in 
the longitudinal direction, and it is rotation-invariant in the transverse direc- 
tion, but x is not boost-invariant in the transverse directions. Hence we have 
generated cylindrically symmetric, longitudinally boost invariant solutions of 
relativistic hydrodynamics. In the longitudinal direction, these solutions are 
homogeneous, boost-invariant and also scale- invariant. Due to this reason, the 
observable rapidity distribution is 



a flat distribution, corresponding to the ultra-relativistic nature of the solution 
in the longitudinal direction (where y = 0.5log[(E + k z )/(E — k z )\ is the 
rapidity of a particle with four-momentum (E, k) and dn/dy is the rapidity 
distribution of particle density). 

A new hydrodynamical solution is assigned to each non-negative function 
V(x), similarly to the cases of the non-relativistic solutions of ref. [13] and 
the 1+1 dimensional relativistic solution of ref. [18]. Note that the solutions 
are valid also for massive particles, the form of the solution is independent of 
the value of the mass m. The form of solutions depends parameterically on k, 
that characterizes the equation of state. 

We have obtained new solutions of the (1+3) dimensional relativistic hydro- 
dynamical equations which describe a self-similar, streaming flow. In the case 
of R = 1 and V(x) = 1 we recover the spherically symmetric scale-invariant 
solution. This means that, in this limiting case, the pressure, the density and 
the temperature profiles depend only on the proper time r. 

In the general case, our solution depends on the scale R (or we can say it 
contains some characteristic scale R) and also on an arbitrary scaling function 



dn 

— = const, 

dy 



(38) 



V(x). 
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4 Summary 



We have found a new family of solutions of 1+3 dimensional relativistic hy- 
drodynamics with cyllindrical symmetry. This family solves the continuity 
equation and the conservation of the energy - momentum tensor of a perfect 
fluid, assuming a simple equation of state, similarly to the case of the recently 
found 1+1 dimensional solutions. The mass of the particles m and the con- 
stant of proportionality between the kinetic energy density and the pressure, 
k, are free parameters of the solution. 

As compared to the well-known case of the scale-invariant solution, we have 
solved one more equation, the continuity equation. We have considered equa- 
tions of state that have two free parameters, the mass m and k = de/dp, 
while the scale-invariant solution is obtained in the m = approximation. In- 
terestingly, our generalizations resulted in additional freedom in the solution. 
The new solutions, similarly to the scale-invariant case, prescribe scaling 3- 
dimensional flow and pressure distribution. However, in our case, the pressure 
is a product of the local number density and the local temperature, hence one 
of these can be chosen in an arbitrary manner. 

The essential result of our paper is that we found a rich family of exact ana- 
lytic solutions of relativistic hydrodynamics that contain both a longitudinal 
Bjorken flow (that is frequently utilized in estimations of observables in high 
energy heavy ion collisions) and a relativistic transverse flow (whose existence 
is evident from the analysis of the single particle spectra at RHIC and SPS 
energies) . 
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